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ABSTRACT 

The radial density profile of dark matter structures has been observed to have an almost universal behaviour in 
numerical simulations, however, the physical reason for this behaviour remains unclear. It has previously been 
shown that if the pseudo phase-space density, p/cr e d , is a beautifully simple power-law in radius, with the "golden 
values" e = 3 and d = r (i.e., the phase-space density is only dependent on the radial component of the velocity 
dispersion), then one can analytically derive the radial variation of the mass profile, dispersion profile etc. That 
would imply, if correct, that we just have to explain why p/a 3 r ~ r~ a , and then we would understand everything 
about equilibrated DM structures. Here we use a set of simulated galaxies and clusters of galaxies to demonstrate 
that there are no such golden values, but that each structure instead has its own set of values. Considering the 
same structure at different redshifts shows no evolution of the phase-space parameters towards fixed points. There 
is also no clear connection between the halo virialized mass and these parameters. This implies that we still do 
not understand the origin of the profiles of dark matter structures. 

Subject headings: - galaxies: halos - dark matter - methods: data analysis - methods: numerical - 



1. INTRODUCTION 

According to numerical simulations of dark matter (DM) struc 
tures, the mass density profile, p(r), changes from something 
with a fairly shallow profile in the central region, 7 = dlnp /dlnr^ 
-1 (or maybe zero), to some thing steeper in the outer region, 

X 



-3 (or maybe steeper) dNavarro et al. 19961: [Moore et al 



1 998 tjDiemand et al^2007l) ( see alsolReed et al l d2003l):ISt oehr 



d2004h ; [Navarro et al.l(l2004l):lGraham et al.ld2006l) : lMerritt et all 
(20061): lAscasibar & Gottloeberl d2008l) ). For the largest struc- 
tures, like galaxy clusters, there appears to be fair agreement 
between numerical pred ictions and observations concerning the 
central steepness dPointecouteau et al.l | 2005b ISand et al~l 20041; 
iBuote & Lewisl2004l:lBroadhurst et alj2005tfvikhlinin et alj2006h . 
however, for smaller structures, like galaxie s or dwarf galaxies , 
observations tend to indicate central cores dSalucci et al.|[2003t 
iGilmore et a"T]|2007t IWilkinson et"aT1l2004l) . Few purely theo- 
retical attempts hav e been made to understand the o r igin of this 
densi ty profile, e.g. Gonzalez-Casado et al. l d2007h : lHenriksenl 
(2007), with varying level of success. 

A completely different approach is to search for simple phe- 
nomenological relations in the numerical simulations, such as 
finding straight lines in some parameter space. The idea is then 
that such phenomenological relations may reduce the complex- 
ity of the Jeans equation, which can then be solved analytically. 

One of the most successful attempts in this direction was 
sparked by the discovery that the pseudo phase-space density 
is approximately a power-law in radius, p/a^ ~ r~ a (Taylor & 
Navarro 2001). The most simple analytical solutions to this 
problem show ed, that the den sity slopes could vary in the range 
from -1 to -3 (Hansen 2004), in excell ent agreement with nu- 
merical results of Nava rro et al.l dl996l) . The analytica l inves - 
tigations were taken to a higher level in lAustin et al.1 J2005), 
where it was demonstrated that there is a characteristic value 
a = 1 .944 when one considers iso tropic structures. Shortly af- 
ter lDehnen & McLa ughlin (2005) used the results of numerical 
simulations dDiemand et al.ll2004alfbT) to show that the "golden 
values" a = 1 .944 and e = 3 indeed provides a very good fit, 
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when one is using the radial velocity dispersion in the pseudo 
phase-space density. Dehnen & McLaughlin (2005, henceforth 
DM05) also solved the Jeans equation under this assumption, 
and demonstrated explicitly that one hereby can derive analyti- 
cally all relevant profiles for the DM structure. Many other au- 
thors have considered similar pseudo phase-space densities, e.g. 



Hansen et alJd2006b:lKnebe & WieBnerid2006h:IStadel et al.ll20 08): 
Knollmann et alJd2008l).lAscasibar & Gottloeberl d2008l):IZait et al.l 
d2008|) : IVan Hese et alJd2008l) : lNavarro et al.ld2008l) : lLapi & Cavaherel 
(2008). 



All this implies, that if we can explain the origin of the very 
simple connection, p/aj ~ r~ a , then we have complete under- 
standing of the DM structures. However, this is under the im- 
plicit assumption that the 3 golden values are indeed the same 
for all structures, namely that a = 1.944, and that the relevant 
quantity to consider is the radial dispersion, erf, with d = r and 
e = 3. 

We will here use the results of recent numerical simulations 
to demonstrate that this is not the case, and that there is no 
simple universal pseudo phase-space density for equilibrated 
DM structures. Given our findings it therefore appears that 
few theoretical approaches that successfully explain the ori- 
gin of the cosmological profiles such as the Barcel ona model 
(Manrique et al. 2003; Gonzalez- Casado et al . 2007) remain. 

2. GENERALIZED PSEUDO PHASE-SPACE DENSITY 

In order to test whether a generalized phase-space density 
exists, we consider the relation 
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(1) 



Here we have defined the general velocity dispersion as dHansenl 
2007; Schmidt 2008) 



e/2 



(2) 



Here (3(r) = 1 - 3r* is me usual velocity anisotropy parameter, 
where a T& & and a Vdn is the radial and tangential component of 
the velocity dispersion respectively. Thus setting D = 0, cor- 
responds to using the radial component of the DM structure 
velocity dispersion in the phase-space density expression. Al- 
lowing D 7/ the phase-space density depends on a velocity 
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dispersion which can be any combination of (T ra d and er tan , e.g. 
D = —i corresponds to using a tot , and D = 1 corresponds to er tan . 

Here we use the rather simple analytical pseudo phase-space 
density, however the actual 6 dimensional phase-space density 
is different from this o ne and is not a p ower-law in radius ac- 
cording to simulations dStadel et al.ll2008h . 

3. NUMERICAL SIMULATIONS 

To test if the values of D, a and e are the same for all struc- 
tures, we used a set of intermediate and high resolution sim- 
ulated DM structures. These structures are all created using 
the PKDGRA V treecode by Joachim Stadel and Thomas Quinn 
dStadelll2001l) . The high resolution simulation 'Via Lactea' in- 
cludes one highly equilibrated structure of mass M; M / = 1 .77 x 
10 12 M p), containing about 84 million particles ( Die mand et al.l 
2007). This structure did not experience any major mergers 
since z = 1 and all the quantities are extracted in spherical bins. 
The rest of the structures are galaxy-size and cluster-size DM 
halos based on either a WMAP 1 year or a WMAP 3 year cos- 
mology. The initial condition s for these structures are generated 
with the GRAFTC2 package (Be rtschingerll2001l) . The starting 
redshifts z, are set to the time when the standard deviation of the 
smallest density fluctuations resolved within the simulation box 
reaches 0.2 (the smallest scale resolved within the initial con- 
ditions is defined as twice the intra-particle distance). All the 
halos were identifie d using a spherical overdensity algorithm 
dMaccio et al 12007T) . The cluster-like halos have been extracted 
from a 63.9 Mpc//i simulation containing 600 3 particles, with 
a mass resolution of m p = 8.98 x 10 7 M q //l The masses of the 
clusters used for this study are 2.1, 1.8, and 1.6 x1O 14 M //j. 
The galaxy-size halos have been obtained by re-simulating ha- 
los found in the previous simulation at high resolution. The 
simulated halos are in the mass range 0.9-2.5 x 10 12 M Q jh and 
have a mass resolution of m p = 4.16 x 10 5 Mq//z. That gives a 
minimum number of particles per halo of about 2.5 x 10 6 . The 
high resolution cluster Chr- W3 has 1 1 million particles within 
its virial radius and a mass of M= 1.81 x 10 14 Mq//z. 

From these numerical simulations we directly calculate all 
the relevant quanteties, such as p(r), o>(r), ag(r), <7</,(r) where 
the (t's are combined to obtain (3(r). These profiles can then be 
compared to the pseudo phase-space density defined in Eq. ([1]). 

4. MONTE CARLO CODE 

I n order to test wheth er the suggested golden values of DM05 
and lAustin et alj (2005) do indeed exist, we wrote a Monte 
Carlo (MC) code to optimize the parameters of the phase-space 
density in Eq. (Q]), for each of the simulated DM structures. 

The MC code is based on the tempera ture annealing principle 
dKirk patrick et al. 1983: iHansenl 120041) . We want to optimize 
the parameter set (D, a, e), so that it makes the LHS and RHS of 
relation (Q~|i converge towards the expected power-law relation. 
In order to do that we search the parameter space and for each 
jump estimate the x 2 value of the relation defined by 

X = £« I "T7"7 \ — • ( 3 ) 

V df 2 (xi) J 

Here xi corresponds to the data input from the simulations, f\ 
and fi corresponds to the LHS and RHS of the relation and dfi 
is the error on /?. Since we are dealing with simulations we 
have no reasonable estimate of the error dfi. Therefore we use 
dfo = 0.054. Choosing different kinds of errors (e.g. -4-, 

■4 and r~ a ) with different magnitudes (0.05, 0.07 and 0. 10) has 
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FIG. 1 . — The optimized D and a values from the MC code plotted against 
each other. The vertical an d horizontal dashed lines indicate the golde n val- 
ues suggested by Austin et al. 12005) and Dehnen & McLaughlin (2005). The 
solid line is a linear fit to the results for the z = structures (circles) and is 
given by a = (0.19±0.02) X D + (1.94±0.02). This indicates that the sug- 
gested golden values are just a consequence of using o> in the phase-space 
density relation. The results for the galaxy-size WMAP1 structures at z = 0.2 
(diamonds) and z = 0.5 (triangles) are over-plotted for comparison, and shows 
a trend very similar to the z = case. No correlation between halo mass (indi- 
cated by different symbol size) and the fitting parameters is detected. 



no significant systematic effect on the final result. Since there 
is the possibility of local minima in the ' Y 2 -lan dscape' we have 
imple mented the metropolis choice in our code (Metropolis et al.l 
1 1 953b . All technical details of this code can be found in Schmidt 

5. RESULTS 

Combining the simulated DM structures with the MC code, 
we are able to estimate the parameters that optimizes the phase- 
space density relation for each structure. If there should be a 
general phase-space density relation, each structure should have 
the same optimized parameters. We see in Fig. [T] that this is 
not the case. In Fig. Q~|we have indicated the suggested golden 
values as horizontal and vertical dashed lines. 

The obtained (roughly) linear relations in Figs.[T]and|2]are 

a = (0.19±0.02)xD + (1.94±0.02) (4) 
e = (0.97 ±0.37) xD + (3. 15 ±0.29). (5) 



Table 1 

Phase-Space Density fitting param eters. The cas e 

MARKED WITH AN * IS THE ONE USED InIHANSEnI ( |2004|) ; 

IAustin et AL.Id2005l ); lDEHNEN & McLaughlin1 ( 12005I) . 



D 



1 

0* 

-1 

_2 

3 



Phase-Space Density 



p/[2a 2 
Pi 



2 \ ^rh + CT, 

ai + ol + a 



e/2 
/2 



2.13 ±0.03 
1.94 ±0.02 
1.75 ±0.03 
1.81 ±0.02 



4. 12 ±0.47 
3.15±0.29 
2.18±0.47 

2.50±0.38 



dark matter structures were not that trivial 
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FIG. 2. — The same figure as Fig.fT]for D and e with a linear fit given by 
e = (0.97 ± 0.37) x D + (3. 15 ± 0.29). 



This shows that a generalized phase-space density relation does 
not exist. Thus our results suggest that the hunt for a physi- 
cal explanation of the (often assumed universal) power-law ap- 
pearance of p/er 3 , is probably a dead end. It seems that this 
expression is nothing more than a possible fitting function with 
the nice property of having the same physical units as the phase 
space density. 

If we force D = (like DM05) we get from the relations © 
and © that a = 1.94 ±0.02 and e = 3.15 ± 0.29. These values 
are in excellent agreement with the results from DM05. If we 
use other values of D, i.e., phase-space densities with combina- 
tions of the different velocity dispersion components we get the 
parameter values listed in tableQ] 

Knowing that the parameters that optimizes the phase-space 
density relation in Eq. (Q]) for simulated structures at redshift 
z = are related, it would be interesting to see whether such 
relations are also present at higher redshifts or not. We use z > 
snapshots for the WMAP1 galaxy-size structures (clusters are 
not significantly relaxed at high redshifts). Running the MC 
code with these DM halos gave results very similar to the ones 
for the z = structures. In general we did not find any indication 
of a redshift dependence for the optimized values, showing that 
there is no special attractor for the values of D, a and e. The 
results at z > (triangles and diamonds) are plotted together 
with the z = ones in Figs. [T][3] The obtained linear relations 
between the optimized parameters are not affected significantly 
by redshift, and therefore our calculations suggest that Eqs. (01 
and Q are valid for all redshifts. 

Furthermore we find no significant correlation between the 
virial mass of the simulated structures and the a, e and D values 
(see Figs.Q][3]l. 

6. CONCLUSIONS 

Using a set of numerically simulated galaxy and cluster sized 
DM structures and analysing them with a Monte Carlo code, 
we show that no generalized pseudo phase-space density rela- 
tion seems to exist in general. We have thus shown, that the 
previously suggested relation pja 3 ~ r~ a does not hold univer- 
sally. The redshift and mass independence of our results show 
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FIG. 3. — The same figure as Fig.[T]for a and e with a linear fit given by 
e = (6.39 ±1. 44) xa+ (-9.14 ±2.65). 



that there is no special attractor for the parameters describing 
the generalized phase-space density. 

Instead we happen to identify a set of seemingly linear rela- 
tions between the parameters D, a and e (describing the gen- 
eralized pseudo phase-space density from Eq. ([1)), which we 
have parametrized in Eqs. (0 and @. 

Thus, given our findings that p/a 3 is nothing but a nice fitting 
formula and not a physical attractor, we are still far from truly 
understanding the density profile of DM structures. 
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